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Electronic and nuclear spins of shallow donors in Silicon are attractive candidates for qubits in 
quantum computer proposals. Shallow donor exchange gates are frequently invoked to preform 
two-qubit operations in such proposals. We study shallow donor electron properties in Si within 
the Kohn-Luttinger envelope function approach, incorporating the full Bloch states of the six band 
edges of Si conduction band, obtained from ab initio calculations within the density-functional and 
pseudopotential frameworks. Inter-valley interference between the conduction-band-edge states of 
Si leads to oscillatory behavior in the charge distribution of one-electron bound states and in the 
exchange coupling in two-electron states. The behavior in the donor electron charge distribution is 
strongly influenced by interference from the plane-wave and periodic parts of the Bloch functions. 
For two donors, oscillations in the exchange coupling calculated within the Heitler-London (HL) 
approach are due to the plane-wave parts of the Bloch functions alone, which are pinned to the 
impurity sites. The robustness of this result is assessed by relaxing the phase pinning to the donor 
sites. We introduce a more general theoretical scheme, the floating-phase HL, from which the pre- 
viously reported donor exchange oscillatory behavior is qualitatively and quantitatively confirmed. 
The floating-phase formalism provides a "handle" on how to theoretically anticipate the occurrence 
of oscillatory behavior in electronic properties associated with electron bound states in more general 
confining potentials, such as in quantum dots. 

PACS numbers: PACS numbers: 71.55.Cn, 03.67.Lx, 85.30.Vw 

I. INTRODUCTION 

Doping in semiconductors has significant technological impact. As transistors and integrated circuits decrease in 
size, the physical properties of the devices are becoming sensitive to the actual configuration of impurities.^ A striking 
example is the proposal of donor-based silicon quantum computer (QC) by Kane,^ in which the monovalent "^^P 
impurities in Si are the fundamental quantum bits (qubits). This intriguing proposal has created considerable recent 
interest in revisiting the donor impurity problem in silicon, particularly in the Sii^^^P system. 

Two-qubit operations for the donor-based Si QC architecture, which are required for a universal QC, involve precise 
control over electron-electron exchange^ and electron- nucleus hyperfine interactions. Such control can presumably 
be achieved by fabrication of donor arrays with accurate positioning and surface gates whose potential can be precisely 
controUediSiSiLS However, we have showni^ that electron exchange in bulk silicon has spatial oscillations on the atomic 
scale due to the valley interference arising from the particular six-fold degeneracy of the bulk Si conduction band. 
These oscillations place heavy burdens on device fabrication and coherent control, because of the very high accuracy 
requirement for placing each donor inside the Si unit cell, and/or for controlling the external gate voltages. 

The potentially severe consequences of these problems for exchange-based Si QC architecture motivated us and 
other researchers to perform further theoretical studies, going beyond some of the simplifying approximations in the 
formalism adopted in Ref. 9, and incorporating perturbation effects due to applied strain^" or gate fieldsii Both 
these studies, performed within the standard Heitler-London (HL) formalism, essentially reconfirm the originally 
reported difficulties regarding the sensitivity of the electron exchange coupling to donor positioning, indicating that 
these may not be completely overcome by applying uniform strain or electric fields. At this point it is clear that 
the extreme sensitivity of the calculated exchange coupling to donor relative position originates from interference 
between the plane-wave parts of the six degenerate Bloch states associated with the Si conduction-band minima. 
Theoretically, this effect is dictated by the HL description of the two-electron singlet and triplet states, defined as 
properly symmetrized combinations of single-particle ground-state functions, where the phases of the Bloch states are 
pinned at each donor site. 
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Our goal in the present study is to assess the robustness of the HL approximation for the two-electron donor-pair 
states. Specifically, we first examine the single donor properties in more detail by including the complete Si conduction 
band Bloch functions. The calculated single donor electron charge density vividly illustrates the rapidly oscillatory 
(and non-commensurate) nature of the donor electron properties. We then relax the phase pinning at donor sites 
and allow small phase shifts in the plane-wave part of the Bloch functions, which could in principle moderate, even 
eliminate, the oscillatory exchange behavior. Within this more general theoretical scheme, which we call the floating- 
phase HL approach, our main conclusion is that, for all practical purposes, phase shifts are energetically unfavorable 
for both singlet and triplet states. The previously adopted HL wavefunctions are thus found to be robust, and the 
oscillatory behavior obtained in Refs. 9 10 11 cannot be taken as an artifact. 

The present paper is organized as follows: In the next section we review the shallow donor problem in Si, fully 
incorporating the details of the Si band structure. We present ab initio results for the bulk Si conduction-band-edge 
charge densities associated with individual Bloch states and single donor states. In Sec. Illll we consider two substitu- 
tional donors in bulk Si, and introduce the floating-phase HL approximation. The two-particle ground state energy 
is compared with that of the standard (pinned-phase) HL states. We present results for the donor electron exchange 
between P donor pairs in Si in situations of practical interest, consistent with the current degree of experimental 
control over donor positioning'"'-^-^ for the fabrication of Si QC. We also indicate how the floating-phase scheme may 
be useful for different Si-based QC architectures. Concluding remarks are presented in Sec. IIVI This work thus 
provides necessary theoretical support and pictures to anticipated experimental studies on qubit exchange coupling 
in a Si matrix. 



II. SINGLE DONOR IN SILICON 



We determine the donor electron ground state using effective mass theory. The bound donor electron Hamiltonian 
for an impurity at site Ro is written as 

^0 = T^sv + 'Hvo ■ (1) 

The first term, Tisv, is the so-called single- valley Kohn-Luttinger Hamiltonian,^-^ which includes the single particle 
kinetic energy, the Si periodic potential, and the screened impurity Coulomb potential 

For shallow donors in Si, we use the static dielectric constant e = 12.1. The second term of Eq. QJ, Hvo, includes 
the inter- valley coupling effects due to the presence of the impurity potential>i^ 

The electron eigenfunctions are written in terms of the six unperturbed Si band edge Bloch states = U|^{r)e^^^'''^■. 

6 6 

V'Ro(r) = E "^M^ - ^o)M^, Ro) = "^M^ - Ro)«M(r)e'''"-("^-^°) . (3) 
/J— 1 /i— 1 

The phases of the plane- wave part of all band edge Bloch states are naturally chosen to be pinned at Rq, and the 
expansion coefficients, also called valley populations, are real. In this way the charge density at the donor site [where 
the donor perturbation potential Eq. (0) is more attractive] is maximum, thus minimizing the energy for ipTia (r) . In 
Eq. |(2Jl, Ffj_(r ~ Rq) are envelope functions centered at Rq, for which we adopt the anisotropic Kohn-Luttinger form 
(e.g., for ^ ~ z, Fj(r) = exp{— [(a;^ + y^)la^ + z'^ /h'^Y^'^} / wiiaFb) . The effective Bohr radii a and h are variational 
parameters chosen to minimize Esv — {'4>'r.q\'Hsv\'4''R.q) , leading to a = 25 A, 6 = 14 A and Esv ^ ~30 meV when 
recently measured effective mass values are used in the minimization.^ 

The Tisv ground state is six-fold degenerate. This degeneracy is lifted by the valley-orbit interactionSfi^ii^ which 
account for intervalley scattering effects and are included here in TLvo ■ Valley-orbit effects are conveniently represented 
by two types of intervalley couplings Ti-vOfi,!^- For valleys at perpendicular directions (e.g., ^ = x, v = z) we take the 
coupling HvOx^z = while for those in opposite directions (e.g., ^ = z, v — —z), 'HvOz,-z = ^^c(l + <5). Of 

course HvOfi^fj. — 0. Taking Ac = 2.16 meV and S = —0.3 correctly reproduces^" the ordering and relative splittings 
of the lowest energy states manifold for P donors in Si: A ground state of Ai symmetry, followed by a triplet of Ti 
symmetry and by a doublet of E symmetry. In unstrained Si, the nondegenerate Ai ground state corresponds to all 
= 1/a/6 in and its binding energy is Eq =< ■0Ro|Ho|V'Ro >= ^sv — (5 + S)Ac ~ —40 meV, to be compared 
to the experimental valuei^ of —45 meV^ Aiming at the ground state of the system, we restrict our discussion to the 
nondegenerate (Ai-symmetry) ground state, thus all = in © in relaxed Si. For strained Si, also considered 

below, the valley populations change according to the degree of strain.— 
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The periodic part of each Bloch function is pinned to the lattice, independent of the donor site. It can be expanded 
over the reciprocal lattice vectors G: 

^^(r)=^cj^e^«"-. (4) 

G 

We determine the coefhcients Cq for the conduction band edge Bloch states in Eq.l^J from ab initio calculations. 
Electron-electron interactions are described by density-functional theory (DFT) within the local-density approximation 
(LDA)ii2iiS We use the exchange-correlation potential parametrized by Perdew and Zunger^° from Ceper ley- Alder 
quantum Monte-Carlo results for the homogeneous electron gag^^. The interactions between valence electrons and ions 
are described by the ab-initio, norm-conserving pseudopotentials of TrouUier-Martins^^, generated by the FHI98PP 
code.^'^ We use 290 plane waves in the expansion of Eq.Q, up to a maximum kinetic energy of 16 Ry. Calculations 
are performed by the ABINIT codej^i The key ingredients of this code are: (i) An efficient Fast Fourier Transform 
algorithms^ for the conversion of wavefunctions between real and reciprocal space; (ii) The use of iterative minimization 
techniques to solve the Kohn-Sham eigenvalue problem, more specifically an adaptation to a fixed potential of the 
band-by-band conjugate gradient method,-^ and a potential-based conjugate-gradient algorithm for the determination 
of the self-consistent potential. Details of the methodology are described in Ref. We find the equilibrium lattice 
constant of Si at a= 5.41 A and the conduction-band minima at 0.844(27r/a) from F, in close agreement with 
experimental results>i& These values are used in the calculations presented below. 

In Figs.^a) and ^b) we present the electronic probability density Px = \(f>x\'^ obtained from the single conduction- 
band-edge Bloch state (j)x = Ua;(r)e*''^ '". Visually, our results indicate that this state is predominantly formed by 
\px) atomic-like orbitals, although some d— character may also be present, in consistency with the higher degree of 
delocalization for these states as compared to the conduction band states at the F point.-® Of course, as for any Bloch 
state, the probability density is periodic in the fee lattice. It is also interesting to note that, among the 290 plane- wave 
states included in our basis, over 90% of the spectral weight in the plane- wave expansion of Ux{r) in Eq.l^J comes 
from five reciprocal lattice vectors: G = (0, 0, 0), ^(—1, ±1, ±1). These give the five smallest values of jG-t-k^:] since 
kx ~ 0.844-^(-|-l, 0, 0). This same criterion for the five most relevant coefficients Cq applies to each of the other five 
k^-vectorsiS£ 

In Figs.njc) and^d) we show the total charge density X]^=i e I'^mP- Fig.QIc) shows the characteristic antibonding 
signature of the conduction band state, which was also found by Richardson and Cohen^® for the conduction-band 
density at the X-point in Si (thus not exactly at the band edge). The conduction band edge state of Si has been 
previously studied by Ivey and Mieher^ Our ab initio results are in good qualitative agreement with this earlier 
empirical pseudopotential study. 

We analyze the effects of Si band structure on donor electron wavefunction within two models for the conduction 
band edge states of Si: (i) (?I)p,Ro = e''''''^'""^"^; (u) 0^,Ro = u^{r)e^^'^'^'^~^°\ Model (i) corresponds to the free- 
electron single-plane-wave-per-valley approximation adopted in previous studiesAiS*^ In model (ii) band structure 
contributions are fully incorporated. Regarding the electron probability density plotted in Fig.^for model (ii), model 
(i) would have given completely uniform distributions, consistent with taking = 1, i.e., Cq — SG,r for all fj,. 

The effects of the conduction band states of Si on the donor wavefunctions and charge density are well established 
experimentally.'^^ Particularly, the single impurity charge density is not only an interesting physical property by itself, 
but also foretells the oscillatory behavior in two-donor properties such as exchange. Figs.|5|(a) and (b) give the single 
electron charge density [^(r)^ along a (GDI) crystal plane for a symmetrized state at the conduction-band edge of 
bulk silicon, ^(r) = (v^)"^ '/'m,Ro within models (i) and (ii) respectively. Frame (a) shows that interference 

from the six plane-wave states included in model (i) leads to a periodic charge pattern consistent with a simple-cubic 
lattice of lattice parameter 27r/fc^ ^ 1.18a, with a periodicity which is clearly different (and incommensurate) from 
the atomic positions in the lattice, since jk^j is incommensurate with the reciprocal lattice. Of course a different 
interference pattern would result if the plane waves were not all pinned at site Rq. Results for model (ii) given in (b) 
show that additional interference from the Bloch functions u^{r), which are periodic in the fee lattice, further reduce 
the periodicity of the charge density. 

Figures |2Ic) and|2Id) give the charge density |'0Ro(i")P for the donor state in Eq. within models (i) and (ii) 
respectively. The impurity site Rq, corresponding to the higher charge density, is at the center of each frame. It is 
interesting (and somewhat counterintuitive) that, except for this central site, regions of high charge concentration 
and atomic sites do not necessarily coincide, because the charge distribution periodicity imposed by the plane-wave 
part of the Bloch functions is 27r/fc^, incommensurate with the lattice period a. 
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FIG. 1: (Color) (a) and (b) Probability density for the single Bloch state px ~ \4>x'^ in two different crystal planes. Notice 
the \px > atomic-like signature, (c) and (d) Total probability density for the six conduction-band minima, showing a more 
symmetric structure. White dots represent Si sites in the diamond structure and the color scheme runs from purple (low 
density) to red (high density) 



III. DONOR PAIR 



A. Exchange coupling within the Heitler-London approach 



Within the HL approximation, the lowest energy singlet and triplet wavefunctions for two electrons bound to a 
donor pair at sites Ha and Rb are written as properly symmetrized combinations of iPra and V'Rb : which are in turn 
defined in Eq. © 

- — =L==[i^R^(ri)?AR^(r2)± VRB(i-i)^R^(r2)], (5) 
^2(1 ± 6^) 

where S is the overlap integral and the upper (lower) sign corresponds to the singlet (triplet) state. The energy 
expectation values for these states are 

Et^{^t\nM)^2Eo + ^^^^ (6) 

where Eq is the isolated impurity binding energy and Hq and Hi are usually referred to as Coulomb and exchange 
integralsiiSii^ The energy difference J = Et — Eg gives the exchange splitting. We have previously derived the 
expression for the donor electron exchange splitting in RefT^, which we reproduce here: 

J(R) =^|a^na,|2j'^,(R)cos(k^-k,) -R, (7) 

where are the valley populations defined in Eq. (|3J), and J7^i/(R) are kernels determined by the envelopes. These 
are slowly varying functions of R (explicit expressions are given in Ref. llOtl . monotonically decaying with distance 
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FIG. 2: (Color) Frames (a) and (b) give the electron probability densities on the (001) plane of bulk Si for the bottom of the 
conduction band eigenstate corresponding to a symmetric combination of the six degenerate Bloch states at the conduction 
band edge, calculated within models (i) and (ii) respectively. Frames (c) and (d) give the corresponding probabilities for the 
ground state of a donor in Si within the envelope function approximation. The white dots give the in-plane atomic sites and 
the color scheme runs from purple (low density) to red (high density) 



since only the exponential envelopes centered at each donor, but not the Bloch functions, contribute to it. Below we 
make a few observations before we attempt to go beyond the HL approximation. 

Equation jT)) does not involve any contribution from the periodic part of the Bloch functions Q [in terms of 
additional oscillatory behavior in J(R) or additional contribution to the magnitude of J(R)], which therefore may 
essentially be taken as (r) = 1 . This fact has been pointed out by Wellard et aljii and is a consequence of the pinning 
of the u^{r) functions to the lattice, independent of the donor site, and of their fast oscillating nature. These authors 
calculated some HL integrals with G different from the T point, which were originally neglected in Refs. 9.10 31, and 
confirmed numerically that all approximations adopted here (and in Ref. Il0|) are excellent. We therefore conclude 
that models (i) and (ii), though giving quite different electron probability densities as illustrated in panels (c) and (d) 
of Fig. 121 effectively lead to the same results for the exchange coupling within the HL approximation. 

Although the exchange coupling given in Eq. Q should be applicable to any relative position vector R, including 
the effect of small perturbations in the donor sites into off-lattice positions^ it has been pointed out by Altarelli and 
co-worker3^2i^ that interstitial donors in Si may acquire a deep-center character, invalidating the envelope-function 
treatment adopted here. We therefore focus our study on substitutional (thus shallow) donors in Si. 

Figure |3| illustrates a case of practical concern involving unintentional donor displacements into nearest-neighbors 
sites, when donor pairs belong to different fee sublatticesi^^ The open squares in Fig.|3Ia) give J(R) for substitutional 
donors along the [100] axis, while the open triangles illustrate the different-sublattice positioning situation, namely 
R = Ro + Snn with Rq along the [100] axis and Snn ranging over the four nearest-neighbors of each Rq {(Inn = 
\Snn\ =a\/3/4 ~ 2.34 A). The lower panel of the figure presents the same data on a logarithmic scale, showing that 
nearest-neighbor displacements lead to an exchange coupling reduction by one order of magnitude when compared to 
J(Ro). 
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FIG. 3: (Color online) Calculated exchange coupling for a donor pair versus interdonor distance in (a) unstrained and (b) 
uniaxially strained (along z) Si. The open squares correspond to substitutional donors placed exactly along the [100] axis, 
the lines give the calculated values for continuously varied interdonor distance along this axis, assuming the envelopes do 
not change. The open triangles give the exchange for a substitutional pair almost along [100], but with one of the donors 
displaced by djvjv ~ 2.3 A into a nearest-neighbor site. The lower frames give the same data in a logarithmic scale. When the 
floating-phase HL approach is adopted, the results change negligibly; the filled symbols on the lower left frame give examples 
of calculated corrections (see text). 



Our previous studies^S show that the extreme sensitivity of J(R) to interdonor positioning is eliminated for on- 
lattice substitutional impurities in uniaxially strained Si (e.g. along the z axis) commensurately grown over Sii_a;Gea; 
alloys if R remains parallel to the interface x-y plane. The strain is accommodated in the Si layer by increasing 
the bond-length components parallel to the interface and decreasing those along z, breaking the cubic symmetry 
of the lattice and lowering the six-fold degeneracy of the conduction band minimum to two-fold. In this case, the 
valley populations in the donor electron ground state wave function Q are determined from a scalar valley strain 
parameter x^ which quantifies the amount of strain. Figure |3|[b) gives J(R) in uniaxially strained (along z direction) 
Si for X — —20 for the same relative positioning of the donor pairs as in Fig.|2Ia). Notice that the exchange coupling 
is enhanced by about a factor of 2 with respect to the relaxed Si host, but the order-of-magnitude reduction in J 
caused by displacements of amplitude d^N into nearest-neighbor sites still persists as Snn is not parallel to the x-y 
plane. 



B. Floating-phase Heitler-London approach 

1. Formalism 

In Refs. iandd as in the standard HL formalism presented above, it is implicitly assumed that the phases e ^'^f '^o 
in Eq. Q remain pinned to the respective donor sites Rq = Ryi and Rs, as we adopt single donor wavefunctions 
to build the two-electron wavefunction. Although phase pinning to the donor substitutional site is required for the 
ground state of an isolated donor {Ai symmetry) in order to minimize single electron energy, this is not the case for 
the lower-symmetry problem of the donor pair. In order to minimize the energy of the two-donor system, here we 
allow the phases to shift by an amount (5R along the direction of the interdonor vector R = R^ — Ra, so that the 
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single-particle wavefunctions in Eq. ^ become 
and 

All terms appearing in Eq. ^ are now functions of 5R, which we take as a variational parameter here, chosen 
independently as SUg and SRt to minimize Es and -Et (since singlet and triplet states are orthogonalized through 
the spin part of the wavefunction) . A similar ansatz, the so-called floating functions approach, was suggested by 
Hurleji^ as an improvement over HL for the H2 molecule, with the atomic orbitals symmetrically shifted towards 
each other. When the amplitude of the shift is taken as a variational parameter, energy reduction thus obtained 
leads to a significantly better agreement with experiment for the hydrogen molecule total energy.— Since the phases 
in Eq.Q are responsible for the oscillatory behavior of the exchange coupling between donor electrons in Si, this 
more general variational treatment might lead to changes in the previously reportedSiiSiii behavior of the two-donor 
exchange splitting J ~ Et ~ Eg. 



2. One- and two-center contributions 

Adopting the floating-phase forms given in Eqs. (jSJ and ^ in the HL expression ((SJ leads to a modified expression 
for the expectation value of the energy in Eq. © for the singlet and triplet states. The term 2Eq on the right hand 
side of Eq. @ gives the single-particle single-center contributions from both (isolated) impurities, which should be 
taken here as Ea + Eb ~ 2Ea- For the present model Hamiltonian, within the floating phase HL approximation, we 
get 

EAiSK) = (^r^IHaI^rJ (10) 

= Esv [2(cos (j)^ -I- cos cjjy -I- cos cjjz)'^ + S{cos 2(j)x -\- cos 20j^ -I- cos 2(j)z) — 3] , 

o 

where 0^ — ~(f>-fi — • SH.. The parameters and 6 are defined in Sec. m where their numerical values are also 
given. As expected, for SR = the above expression leads to Ea{0) = Eq = Esv ^ (5 + S)Ac', while SR 7^ leads to 

Ea^SR) ^Eo + ^>Eo. (11) 

The correction Aq is positive definite by construction, since the one-particle functions in standard HL (SR = 0) are 
taken as the ground-state wavefunction of the isolated impurity problem. 
The expectation value for the energy of the donor pair is given by 

E't{R, SR) = £;|(R, 0) + Ao((5R) + A'^{R, SR) , (12) 

where i?((R, 0) is the pinned- phase result from the regular HL calculation. The first correction term, Ao((5R), is the 
energy shift due to the single-particle single-center contributions, derived above, and A((R, JR) are the singlet and 
triplet state corrections coming from the two-center contributions Hq, Hi and S. The latter are integrals involving 
the electronic wavefunctions © and ®, and are calculated here as described in Ref. 0, with the proper phase shifts 
included in the plane- wave part of the Bloch functions. 

In Fig. ^ we give the calculated values of the individual corrections Ao(<5R) and A((R, 5R) for a geometry where 
the impurities are 16 lattice constants apart {'--^ 87 A), with R along the [100] crystal direction. The energy correction 
Aq raises sharply for nonzero SR, and is of course independent of the relative position vector R, while the energy 
variations Af oscillate and decrease with increasing relative distance R, and may be positive or negative according to 
SR (for SR ~ in the case illustrated in Fig. ^ As decreases for negative shifts SR, while At decreases for positive 
shifts). Since Aq is always positive, independent of R and very sensitive to SR, we conclude that the effect of phase 
shifts aiming at minimizing two-donor energy is negligible and may be safely ignored for R^ a,b, where a and b are 
the donor effective Bohr radii. For example, minimization of the total energy in Eq. H12|l for the particular geometry 
considered in Fig. pleads to SRg = —7 mA, with the singlet energy decrease of 270 neV, and SRt = -1-7 mA, with the 
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FIG. 4: (Color online) Calculated corrections to the total energy for a P donor pair in Si. The donors are 87 A apart, 
along the [100] direction. The parameter SR is the amplitude of the individual phase shifts from the donor sites, ±(SR, along 
the interdonor line. The solid line gives Ao, the single-center contribution, while the dotted (dashed) line gives As (At), the 
two-center singlet and triplet contributions, respectively. 



triplet energy decrease of 6 neV. This results in an increase in J by (264) neV, given by the solid square in the lower 
left hand side frame of Fig. 13 The floating phases variational scheme leads to a reduction in both singlet and triplet 
states energy, therefore the net variation in J is positive (negative) if the triplet energy reduction is smaller (larger) 
than the singlet. The solid triangle in Fig. |3| corresponds to a case of negative variation, obtained when one of the 
donors in the above geometry is displaced into a nearest-neighbor site. Note that the corrections are more than three 
orders of magnitude smaller than the calculated J assuming 6Rs — SRt = 0. In other words, for all practical purposes 
the fixed-phase standard HL approximation is entirely adequate for the range of interdonor distances of interest for 
QC applications. 

This conclusion is not in contradiction with Hurley's result for the hydrogen molecule, '^^ where significant energy 
reduction is obtained around the equilibrium nuclear separation, R ^ 1.5ao (ao = 0.53 A is the free hydrogen atom 
Bohr radius). For R of the order of the Bohr radius , Ag becomes comparable to Aq, resulting in an improved 
variational estimate for the ground state energy of the H2 molecule when small shifts are allowed in the single-particle 
hydrogenic orbitals. 

Since the current calculation has taken into account the full bandstructure of the host Si material, and modifying the 
standard HL approximation has proved to have minimal effect on the results, further improvement in a perfect crystal 
environment (that is, relaxed bulk Si) can only be achieved by including the higher energy orbitals4i^ However, 
we do not anticipate significant moderation of the fast oscillatory behavior of exchange coupling as all the orbitals 
share the same conduction band valleys, though quantitative shifts might be expected in a larger scale molecular 
orbital calculation. In the present study, we keep the two donors relatively far apart so that the HL approximation is 
applicable. This is also the situation of interest to practical QC fabrication considerations (which requires the donors 
to be at least 100 A apart). 

Another improvement over our current calculation may come from including the effect of lattice distortions. The 
Coulomb interaction between the additional protons on the lattice sites and the two electrons for the donor pair 
creates a strain field on the underlying crystal lattice. Such a field affects the electronic structure in the same way 
as the uniaxial strain discussed above, though it is along the inter-donor axis. Since the inter-donor separation in 
the present situation is much larger than the effective Bohr radius (~ 30 A), the inter-donor interaction is strongly 
screened, therefore lowering the strength of the strain. Furthermore, if a uniaxial strain is already applied along 
the z direction, so that the donor ground state only consists of the z and —z valleySfiS the additional strain due to 
the presence of another donor (e.g. along x direction) will not further reduce the number of valleys involved — the 
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nondegenerate ground state will still consist of an equally weighted superposition of these two valleys (instead of 
just one of them), so that oscillations in exchange due to valley interference cannot be removedflS, Nevertheless, a 
quantitatively analysis is needed to assess the significance of this effect. 

3. Coupled quantum dots 

Shallow donors in semiconductors may be viewed as the simplest, naturally occurring quantum dot. Compared 
to the gated quantum dots, a relevant difference is the presence of a well defined and sharp pinning center at the 
substitutional donor site. Previous proposals of quantum dots as quantum gates in a Si or Ge matrijsSSiffi were based 
on estimates for the exchange coupling within an empty envelope function description. It is clear that, for these 
materials, the plane-wave parts of the Bloch functions may also have an important effect in the exchange coupling. 

The floating-phase HL approach should be applicable to coupled quantum dots, leading to an expression equivalent 
to Eq. (|12|l . The absence of a sharp pinning center associated with each quantum dot implies that Aq is not as 
sensitive to the phase shifts in a floating-phase variational scheme as obtained here for the donor case. It is possible 
that variations in the two-center contributions Af dominate energetically and determine the singlet and tripled ground 
state energies, whose difference should give a reliable estimate for J. 

Of course the valley-orbit effects described by Hyo in SeclHl which are quantitatively well established for P donors 
in Si, would have to be estimated for the quantum dot confining potential, including other perturbations which break 
the translational symmetry of the host potential, such as the presence of nearby interfaces and strain^ As in the 
present case, Hvo should lift the six-fold degeneracy of the isolated quantum dot ground state. An investigation of 
valley-orbit effects in Si quantum wells was performed recently by Boykin el al.^'^ 

A similar scheme may also be useful for spin cluster qubita^ embedded in Si or Ge, where exchange gates are also 
invoked for inter-cluster interactions. Demonstration that the exchange oscillatory behavior is circumvented for spin 
clusters would further require the formalism to be generalized to include multielectron states^ in each cluster, as was 
explored in Ref. H^- 

IV. CONCLUDING REMARKS 

We have included and assessed full band structure effects in the single donor wavefunctions and charge distributions 
in Si. We find interesting oscillatory patterns resulting from interference between the different plane wave components 
of the Bloch functions. Regarding the well-separated donor pair problem, we introduced a generalized scheme — the 
floating-phase HL approach, which reconfirmed the reliability of standard HL for this range of donor separations. 

One perceived advantage^ of Si-based spin quantum computation (over, for example, the corresponding GaAs 
quantum dot based quantum computation) is the universal nature of each qubit in Si, i.e. the fact that the P 
donor electronic state in Si is always exactly the same, making each qubit identical (without any need for additional 
characterization of individual qubits which will surely be needed for GaAs quantum dot quantum computers since 
electrostatically confined electronic spin states in GaAs quantum dots would obviously have a fair amount of qubit 
to qubit variations as no two quantum dots can really be identical). Our finding of exchange oscillations in Si donor 
states demonstrates that this perceived advantage of Si comes with a price, where the exchange coupling between 
qubits may vary depending on the precise positioning of the P atoms within the Si unit cells. We believe that, in 
spite of this problem, the QC scheme with donors in Si still has its appeal in terms of uniform qubits. Obviously 
some characterization of the exchange coupling in Si becomes necessary in view of the oscillatory exchange behavior. 
We have discussed elsewhere^ how some precise local information about donor state exchange coupling in Si can be 
obtained by using the powerful tool of the micro- Raman scattering spectroscopy. In addition, various band engineering 
procedures, "'^^ using strain effects and/or Si-Ge quantum dots, could be utilized to reduce the exchange oscillation 
effects, although its complete elimination may not be easy. 

From the perspective of current QC fabrication efforts, ^ 1 nm accuracy in single P atom positioning has been 
recently demonstrated)^ representing a major step towards the goal of obtaining a regular donor array embedded in 
Si. As expected,^ electronic calculations^ have confirmed that this degree of control is entirely compatible with the 
operations involving the so-called A-gates in the Kane qubit architecturei^ On the other hand, the present calculations 
have confirmed that deviations in the relative positioning of donor pairs with respect to perfectly aligned substitutional 
sites along [100] lead to order-of-magnitude changes in the exchange coupling. Severe limitations in controlling J would 
come from "hops" into different fee sublattices, in particular among nearest-neighbor substitutional sites. Therefore, 
precisely controlling exchange gates in Si remains an open challenge. 
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